Optical Resonances in Reflectivity near Crystal Modes with Spatial Dispersion 
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We study the effect of spatial dispersion of crystal modes on optical properties such as the re- 
flectivity R. As an example for isotropic media, we investigate the simplest model for phonons in 
ionic crystals and compare with previous results for highly anisotropic plasmons, which are now 
understood from a more general point of view. As a consequence of the wave vector dependence of 
the dielectric function small changes in the lineshape are predicted. Beyond that, if the frequency 
of minimal R is near a pole of the dispersionless dielectric function, the relative amplitude of dips 
in R with normal and anomalous dispersion differ significantly, if dissipation and disorder are low. 



Usually, in crystal optics the dependence of the di- 
electric function on the wave vector k is neglected be- 
cause the wave length A is much larger than interatomic 
length scales s and effects of spatial dispersion are ex- 
pected to be of the order of s/A -C 1, see Refs. and 
references therein. Recently, it was shown for the re- 
flectivity near the Josephson plasma resonance (JPR) in 
highly anisotropic superconductors^, which is an inter- 
layer charge oscillation of Cooper pairs, that spatial dis- 
persion provides a novel way to stop light and can have 
effects of order unity on the relative amplitude of plasma 
modes with different dispersion. This system is excep- 
tional with respect to the symmetry (high anisotropy) , 
the fact that light is screened away from the JPR, and 
the type of the excitation (plasmon). Here we discuss 
the possibility that such effects appear under more gen- 
eral conditions. As a generic example of a spatially dis- 
persive resonance we choose the simplest model of one 
optical phonon band in an ionic crystal^ (also describing 
certain excitons 1 ) and study the optical properties of the 
isotropic phonon-polariton, which is the lattice vibration 
coupled to light. Thereby, we determine under what cir- 
cumstances similar effects can be expected in a large class 
of systems and clarify the origin of previous results for 
the highly anisotropic plasmon obtained in Ref. 0. We 
show that the effect of spatial dispersion is of order unity, 
if the resonance in reflectivity is in a frequency region 
of size Alu/ujq = s/A near the pole u>q of the dispersion- 
less dielectric function e(w); the amplitudes in reflectivity 
of phonon bands with normal and anomalous dispersion 
then differ significantly. 

Weak disorder and small dissipation turn out to be es- 
sential for the observation of significant effects due to spa- 
tial dispersion, which is the case in the plasmon system 
of the JPR below the superconducting gap. In insula- 
tors, e.g. in ionic crystals with perfect lattice symmetry 
and weak disorder, such favorable conditions might be 
realized as well, while in metals similar effects would be 
strongly overdamped. 

In the following, we mainly consider an electromag- 
netic wave of frequency uj in normal incidence, with a 
wave vector k = (0, 0, k z ) at the surface z = of a semi- 
infinite dielectric medium. The macroscopic Maxwell 
equations describing the electromagnetic fields inside the 



medium expressed in Fourier components take the form 
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where E x and B y are the components of the electric 
and the magnetic fields, respectively. For an isotropic 
medium the normal modes are determined by 



(2) 



For definiteness, we consider the spatially dispersive di- 
electric constant e(u),q), q = sk z , q G [0, 2ir] (s lattice 
constant) , of an ionic crystal, which can be obtained by 
solving the microscopic equation of motion for the ionsi 
(using the force constants of the bulk) , 
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where is the background dielectric constant, S 



S (l + aq 2 ), a ~ KT 3 , S - 
lu 2 = u 2 (l+dq 2 ), wo/(2tt) 



1, is the oscillator strength, 
- 1 THz, d ~ 0.1, is the dis- 



persion of the crystal mode, and p is the damping due 
to disorder (p/ljq ~ 1CP 2 for NaCl, see Ref. Since 
a/d -C 1, the wave vector dependence of the oscillator 
strength can be neglected. The sign of d determines the 
shape of the dispersion and both d > (normal disper- 
sion) and d < (anomalous dispersion) naturally appear 
in ionic crystals due to the backfolding of the modes into 
the first Brillouin zone (e.g. realized in fee structured 
NaCl for the directions k|| [100] and k||[lll], see Ref.0). 
The dielectric function in Eq. J2J is also a model for cer- 
tain excitons where the damping p/loq ~ 10~ 4 is much 
smaller, see Ref. Q. 

The formula analogous to Eq. J2J) for an anisotropic 
plasmon such as the JPR reads 



+ 



/i . 



e c (w, k z ) e a (w) 



(4) 



the 



spatially 

,2/ 



dispersive 

,2\ 



dielectric 
,2/ 



function 

, ,2 



where 

e c {u,q) = e o0 (l - uj 2 c {q) / uj' 2 ) , q = sk Zl uj 2 {q) = uj 2 q (1 
dq 2 ), accounts for the c axis plasmon (excitation paral- 
lel to the z axis). This mode is only excited for light in 
oblique incidence k = (k x ,0,k z ), where k x = ujsa\6/c 
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FIG. 1: Polariton dispersion relation co(\q\) for a transverse 
optical (TO) phonon. Multiple modes qi,q2 are excited by 
incident light of frequency uj according to the solutions of 
Eq. @. The appearance of an extremal frequency cu e , where 
the group velocity v gz — dui/dk z vanishes and q changes sign, 
is characteristic for the mixing of normally dispersive light 
and a crystal mode with anomalous dispersion. 
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FIG. 2: A pole u>o in k 2 (uj) in the case without spatial dis- 
persion (dashed lines) is regularized when spatial dispersion 
is taken into account (solid), see also the inverse uj(k 2 ) in 
Fig-CHa). Depending on the sign of d oc dk 2 (i>j)/dw, an ex- 
tremal point uj e , where the group velocity v gz — du>/dk z van- 
ishes, occurs in case of anomalous phonon dispersion (a) and 



an extremal point an, where k 2 (uji) 
case of normal phonon dispersion (b). 



-k 2 (uJi), appears in 



is conserved due to translational invariance parallel to 
the surface (9 is the angle of incidence). The dielectric 
constant parallel to the surface e a (tu) = e a o(l — w a/ w2 ) 
is non-dispersive. The resonance frequency u a for the 
oscillation parallel to the surface is far above the c axis 
plasmon frequency, uj a 3> w C 0i an d for frequencies near 
u! c o it follows that e a (u) ~ —Lo^juJ 1 <C —1. 

Focusing on the phonon case first, the dispersion rela- 
tion ui(q) (see Fig. ^| follows from Eq. © with e(q,cu) 
from Eq. @ and shows that at a given frequency w 
two modes with (generally complex) refraction indices 
rip = cq p /suj, p — 1,2, can be excited due to spatial 
dispersion. This corresponds to the regularization (cf. 
Fig. I2J) of the pole in k 2 (ui) oc e(u), cf. Eq. ©, which 
exists in the dispersionless dielectric function e(w) at the 
phonon resonance loq in the absence of damping (p — 0) . 
Thereby the appearance of an extremal frequency u> e , 
where the group velocity v gz = duj/dk z for propagation 



in the z direction vanishes, is characteristic for the mix- 
ing of the normal dispersive light mode and an anomalous 
crystal mode in Fig.^a) (cf. Fig.QJ, while for equal dis- 
persion of light and crystal modes a special frequency u>i 
appears, where qf(uJi) — ~q%(u)i), see Fig. Near 
u> ei i the amplitude and lineshape in the reflectivity can 
be sensitive to the type of dispersion, see below. 

In a semi-infinite crystal, the causality principle re- 
quires that v gz is positive, thus only two of the four so- 
lutions of Eq. J2J for n p are physical, e.g., for an anoma- 
lous dispersion the real part of n p is negative. When 
^ l n 2|j only the small-g mode ri\ is light like and 
expected to affect optical properties, while for |ni| « |ri2| 
both modes have to be considered. Note that for two real 
modes near u e , where v gz = (Fig.Q, causality demands 
ni = -n 2 . 

The excitation of multiple modes poses the problem 
that the usual Maxwell boundary conditions, requiring 
continuity of the electric and magnetic fields parallel 
to the surface, are insufficient to calculate the relative 
amplitudes of these modes. In the past, this problem 
was addressed in a phenomenological way by introducing 
so called additional boundary conditions (ABC) for the 
macroscopic polarization P. We use the general ABC 1 



P(*) + 1 (BP/dz) = (z-»0), 



(5) 



where / is a phenomenological decay length and P = 
x(w,k)E, with x^k) = (c(ui, k) — f m ) /An for the 
isotropic medium. These ABC have been justified for 
the anisotropic plasmon in a microscopic model by sub- 
tracting the true equation of motion for the surface layer 
from the hypothetical bulk equation for this layer, see 
Ref. y. The analogous analysis for the phonon allows 
to identify the parameter I in Eq. (J5J with the nearest 
neighbor separation s/2 in fee cubic crystals (e.g. NaCl). 

Given the boundary conditions on E and B, we de- 
termine the reflectivity R at the boundary z = to the 
vacuum, which is the ratio of the reflected and the inci- 
dent energy flux density, 

1-k E x (z = 0) 



R = 



1 



cos 6 B y (z = 0)' 



(6) 



In the dispersionless case the reflectivity R shows a dip 
with a minimal value i? m i n = at e(tu) = k = 1, while 
the frequency range of total reflection, R = 1, between 
the pole and the zero of e is bounded by sharp edges 
in R(oj), see Fig. [21 In the following we will show how 
spatial dispersion can reduce the amplitude 1 — i? m i n of 
the dip in R near k — 1, and modify the lineshape of R 
close to the region of total reflection (cf. Fig. 01 . 

For a spatially dispersive medium in the simplest case 
of normal incident light, where two transverse modes 
q p = (0, 0, g p ), p = 1,2, are excited, we expand the elec- 
tric field E x — Yl p =i Ex(q P )e tq " z ^ s and the polarization 
field P x = El=iEMx(to,q p yi» z / s . Using Eqs. ©, 
(0, and © we derive 
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FIG. 3: (a) Polariton dispersion relation co(q 2 ) for different crystal dispersion d, where q 2 > (q 2 < 0) corresponds to 
propagating (decaying) modes, (b) Dielectric function e(uj) in the absence of spatial dispersion, (c) Reflectivity R(ui) of the 
phonon band. In the region Auj/uiq — (ujq — co e )/ujo ~ s/X near u>o the lineshape is affected by spatial dispersion, but the 
amplitude I — i? m i n of the dip away from this frequency region is not modified. 



in the isotropic system near the extremal frequencies uj ej i 
(ABC for P x ) with the parameter £ = Ico/c = s/2X ~ 
1CP 5 . Similarly, for the anisotropic plasmon with the 
incident field E in the ccz-plane^ (ABC for P z ), we get 



1 



nin 2 



e a cos 9 n\ + ri2 — i£nin2 



(8) 



The above expressions for k enter the reflectivity R, see 
Eq. JfJJl. The difference in Eqs. Q, © arises from the 
anisotropy e a =/= e c in the plasmon system, which can 
suppress the value of n for a strong anisotropy |e | 3> 1. 
Away from w e ,j, where |m| -C |ri2 1 , Eq. @ reduces to 
Ki SO = n\ and Eq. JHJ to Kaniso = ni/(e a cos 9), i.e. the 
usual Fresnel theory is recovered and the smaller (light 
like) refraction index determines the optical properties. 

For d < 0, in a frequency interval Alu/loq ~ y/\d\s/\ ~ 
£ below w e , where ni w —77,2 and the group velocity 
i> ff2 vanishes, the standard Fresnel theory is invalid and 
small length scales, which enter the reflectivity R via 
the parameter £, are relevant. Note that at u e , where 
the (positive) solutions qi of Eq. (0) merge (cf. Fig. 
and ni = —rt2 due to causality, Ki so in Eq. Q is purely 
imaginary and thus R = 1, see Eq. (jSJ. Slightly away 
from a> e , the amplitudes 1 — i? m i n of the dip for isotropic 
phonons and anisotropic plasmons equal to (l+£ 2 ) -1 ~ 1 
and 2/ ((1 + e^£ 2 cos 2 9) 1/2 + l) < 1, respectively. As 
e Q £ ~ 0(1), 1 — i? m i n is significantly suppressed in the 
anisotropic system, while £ <C 1 is negligible for the 
isotropic one. In addition, the band edge near uo e is 
slightly shifted to lower frequencies in comparison with 
the dispersionless theory, (cu — u e )/uu Q ~ £ (Figs.|3©. 

For d > near u>i and above ojq, one solution k 2 of 
Eq. (0) is positive, Fig. |3fa), and thus creates a tail, 
where R < 1 instead of R — 1, see Fig. 0] The difference 
in reflectivity between normal and anomalous phonon 
dispersion for damping p/luq ~ £ is ~ 1% near the band 
edge wo, see Fig. |3Ja). 

These effects on the phonon lineshape in the reflectiv- 
ity R are small and, in particular, the amplitude 1 — i? m i n 
of the dip above u>o is hardly affected, see Fig. |3^c). In 
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FIG. 4: Lineshapes of the reflectivity R near the upper band 
edge of the lower polariton branch for phonons with differ- 
ent dispersion d in case of very low damping p/uo ~ 10~ 6 
(too = 2, So = 3). The shift of the band edge for anomalous 
phonon dispersion is of the order Alu/ljo ~ C ~ 10~ 5 (dis- 
tance between the solid and the dashed lines). In the region 
close to the phonon resonance u>o, where both modes 711,712 
are important, the lineshape depends on the sign of d. 
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FIG. 5: Reflectivity R of the plasma band in the highly 
anisotropic system of Ref. .'! (schematic) . The minimum Rmi n 
near k — 1 is in a frequency region Alu/uo = {u> — tj e ) ~ £ 
close to the bandedge due to the high anisotropy \e a \ 3> 1 and 
spatial dispersion can affect the amplitude of the dip 1 — J? m i n . 
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FIG. 6: Dielectric function e(ui) with usual background di- 
electric constant ~ 1 and in the presence of a back- 
ground mode too = e P oie(<^)- In the latter case the minimum 
Rmin at k — \fe = 1 is shifted to the critical region of size 
Aw/wo = (w e — ujo)/uq ~ £ around the pole u>o (cf. Fig. 
where spatial dispersion can affect the optical properties. 



comparison with the anisotropic plasmon, where a strong 
dependence of the dip amplitude on spatial dispersion 
was reported^, we seek to identify the crucial difference 
to the present system and under what circumstances such 
an effect might appear in the phonon system. 

Firstly, we conclude from Eq. that in an anisotropic 



system (e a ^ e c ) a pole in q 



(e a /e c )sin 



like in Fig.^is created in oblique incidence (9 ^ 0) at the 
zero of e c in a similar way as it appears in the isotropic 
case at the pole of e(uj). More importantly, the strong 
anisotropy \e a \ ^> 1 (~ 10 4 for JPR) of the order ~ l/£ 
suppresses the value of /t a niso m Eq. ©, such that the 
dip in R near k = 1 is close to the bandedge and thus in 
the frequency range Alo/luq ~ £ near the special points 
uj e .i (Fig. where spatial dispersion is relevant (Fig.^J. 

Guided by this observation, we find significant ef- 
fects on the dip amplitude 1 — R m i n for the isotropic 
phonon, if the minimum in R is within the frequency 
interval Aw/wo ~ s/X near the pole of e(u>). This can 
be achieved (a) by a second background dispersionless 
phonon, e w -> e po i c (w) = f m + S b u>l/(u> 2 - u\ - ipbu), 
which is sufficiently close to u>o, (loq — bJb) / u$ ^ s / X ^ 
i.e. |e po ie(wo)| ~ 1/^ ^> 1 and creates an additional 
dip in R, see Fig. [5J or similarly (b) in extreme oblique 
incidence (9 ~ it/ 2 — £}^ 2 )-, where the dispersionless 
K = n/ecos6> = (e - sin 2 6>) (1/2) / e cos6> is I at uj + O{£). 



This results in significant differences in i? m i n for normal 
and anomalous dispersive phonon bands, if the dissipa- 
tion is small, see Fig. [7] Physically the larger suppression 
of 1 — i? m i n for dispersion d > I is related to the fact that 
near u>i only a single propagating mode is excited, while 
the other one decays into the crystal. This effect is of 
order unity and the dip of width Aw /u>o ~ £ in R can be 
resolved in optical phonon spectra™. 

In conclusion, for isotropic phonons small effects for 
the lineshape in reflectivity due to spatial dispersion are 
predicted (Fig.^l. For low dissipation the dip amplitude 
can depend significantly on the sign of the dispersion, if 
the frequency of minimal reflection is in a small distance 
1 
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FIG. 7: Reflection amplitude R for a phonon with different 
dispersion d and small dissipation (a) in the presence of a 
dispersionless background mode eoo — * e po ic (St = 1, Coo = 3, 
p/u)o ~ Pb/uo ~ 10 -7 <C near the band edge at cut = 
0.9995tJo, or (b) for almost parallel incidence (9 w n/2 — ^ 1 ^ 2 , 
p/uia ~ £). In (b) for dispersion d < the minimal reflectivity 
Rmin — 0, but the minimum is slightly shifted compared to 
the dispersionless model (d ~ 0), while for d > the dip 
amplitude 1 — -Rmin is strongly suppressed. 



Aw/wo ~ s/X to the dispersive resonance (Fig. 01 . This 
can be realized in strongly anisotropic systems such as 
the JPR 3 or in the isotropic case due to a second almost 
degenerate resonance or for nearly parallel incidence. 
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